Localization of two interacting electrons in quantum dot arrays driven by an ac-field 
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We investigate the dynamics of two interacting electrons moving in a one-dimensional array of 
quantum dots under the influence of an ac-field. We show that the system exhibits two distinct 
regimes of behavior, depending on the ratio of the strength of the driving field to the inter-electron 
Coulomb repulsion. When the ac-field dominates, an effect termed coherent destruction of tunneling 
occurs at certain frequencies, in which transport along the array is suppressed. In the other, weak- 
driving, regime we find the surprising result that the two electrons can bind into a single composite 
particle - despite the strong Coulomb repulsion between them - which can then be controlled by the 
ac-field in an analogous way. We show how calculation of the Floquet quasienergies of the system 
explains these results, and thus how ac-fields can be used to control the localization of interacting 
electron systems. 

PACS numbers: 73.23.-b, 72.15.Rn, 03.67.Mn, 42. 50. Hz 



I. INTRODUCTION 

Recent years have seen considerable progress in the 
theoretical and experimental investigation of electron dy- 
namics in mesoscopic systems. Preservation of the quan- 
tum coherence of the electronic states can give rise to 
unusual transport phenomena, one of the most notable 
of which is termed "coherent destruction of tunneling" 1,2 
(CDT) , or alternatively "dynamical localization" . This 
occurs when the tunneling of a particle between potential 
minima is suppressed, thus localizing the particle, under 
the influence of an oscillatory electric field. The degree of 
suppression depends sensitively on the parameters of the 
driving field, thereby giving the attractive possibility of 
coherently manipulating electronic states in semiconduc- 
tor nanostructures by means of laser pulses or oscillating 
gate potentials. 

The presence of CDT was first recognized in quantum 
well systems^ and superlattice o 3 : 4 : 5 , in which the effects 
of inter-electron interactions are negligible due to the 
high concentration of carriers. More recent work 6,7,8,9 
has focused on quantum dot (QD) systems, which give 
the ability to study a small number of electrons in a 
highly-controllable environment. In such nanostructures, 
however, the low electron-density means that correla- 
tions produced by the Coulomb interaction cannot be 
neglected, and indeed can have extremely significant 
effects 10,11 . Numerical investigations of a system of two 
electrons in a double QD have revealed novel forms of 
CDT, resulting from the inter-electron Coulomb interac- 
tion, which have been successfully analyzed by means of 
the Floquet approaches. Here we extend and generalize 
these results by investigating the case of two electrons 
confined to a one-dimensional array of QDs, and make 
a detailed study of the interplay between the driving 
field and the Coulomb interaction on coherent transport. 
Although interesting in its own right, such understand- 



ing is also of practical importance to the field of quan- 
tum computation^, as this type of structure provides 
a promising method of implementing scalable arrays of 
quantum bitsi 3 -. 

In this investigation we first demonstrate that the 
Coulomb interaction produces two regimes of behavior, 
namely when the driving is stronger or weaker than the 
Coulomb interaction. This was seen previously in the 
double QD systemic, and we show that it is a generic 
feature, common to all sizes of array. In the strong 
driving regime CDT occurs when the field frequency is 
in resonance with the Coulomb energy at certain field 
strengths which can be predicted by Floquet theorjiiS. 
In the weak-driving regime, however, no resonance ef- 
fects arise, and initial states in which a QD is doubly- 
occupied show a surprising behavior 9, isii£ in that despite 
the strong Coulomb repulsion between them, the two 
electrons can bind together to form a single composite 
particle of charge 2e. The localization and tunneling dy- 
namics of this two-electron state can then be controlled 
by suitably choosing the parameters of the ac-field, in a 
similar way to the control of a single particle. 



II. MODEL AND METHODS 

We consider a system of two-electrons confined to a 
one-dimensional array of QDs, with a time-dependent 
electric field F(t) applied parallel to the array. Such an 
array can be produced experimentally by, for example, 
gating a two-dimensional electron gas confined in a semi- 
conductor heterostructur o 13 : 16 . We model the QD array 
as a single-band tight-binding model of Hubbard type. 
For the double QD case, explicit calculation has shown 
that such an effective model is able to reproduce well the 
behavior of a more realistic modelifi, and accordingly for 
an array of N QDs we use the effective lattice Hamilto- 
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Here U is the Hubbard {/-term giving the energy cost 
for double-occupation of a QD, and t is the hopping pa- 
rameter between adjacent QDs. In experiment t can be 
set by regulating the height of the inter-dot tunneling 
barriers^, and U by adjusting the size of the QDs, since 
U ~ 1/L for a QD of size L. We shall measure all ener- 
gies in units of t, and set h = 1. The operators Cj CT /cj CT 
are the electronic annihilation/creation operators, rij a is 



the usual number operator, and rij 



The 



inter-dot spacing is denoted by a, and for simplicity we 
shall describe the time-dependent electric field in terms 
of the potential difference between neighboring sites, 
E(t) = eaF(t), which we parametrize as: 



E(t) = Esinut, T = 2ir/u. 



(2) 



Experimental investigations of semiconductor QDs 
have revealed that typically the spin relaxation time 
within these structures is very longii and accordingly 
we do not include any spin-flip terms in the Hamiltonian 
QJ. As a result the singlet and triplet sub-spaces of the 
model are completely decoupled. In this work we con- 
sider just the singlet sub-space, as in the triplet sub-space 
the Pauli principle forbids double-occupation of a QD, 
and the Coulomb interaction described by the Hubbard- 
U term is consequently irrelevant. Suitable basis states 
for the spatial component of the two-electron wavefunc- 
tion can be obtained by taking symmetric combinations 
of single-particle states defined on the lattice sites m and 



I , 0m(n.)</>n(r 2 ) + <Mri)0 m (r 2 ) , f9 \ 
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where r±, r 2 are the coordinates of the two electrons. It 
can be clearly seen that the basis consists of N(N — l)/2 
entangled states (Eq[3J) in which the two electrons occupy 
different QDs, with the remaining N factorisable states 
(EqBJ being those in which a QD is doubly-occupied. 
For strong values of the Coulomb interaction the ener- 
gies of these two types of states will fall into two bands, 
separated approximately by the Hubbard gap U. 

Our numerical investigation consists of initializing the 
system in a given state, and propagating it forward in 
time under the action of the time-dependent Hamiltonian 
Q over many (typically between ten and one hundred) 
cycles of the driving field. This can be done efficiently 
by a fourth-order Runge-Kutta method. To study its 
time evolution quantitatively, we measure the projection 
of the two-electron wavefunction onto the basis states 



together with the distribution of electronic charge within 
the array, and at all times we check the normalization 
of the wavefunction to verify that unitarity is adequately 
preserved by the numerical procedure. 

As the Hamiltonian (|TJ is periodic in time, H(t) = 
H(t+T), the Floquet theorem allows us to write solutions 
of the time-dependent Schrodinger equation as = 
exp[— iejt\\(j)j{t)), where \<f>j(t)) are termed Floquet states 
and Cj are the Floquet quasienergies. The Floquet states 
have the same periodicity as the Hamiltonian, and are 
eigenvectors of the unitary time-evolution operator for 
one period of the driving: 



U(t + T,t ) = T exp 



r*o+T 
-i I H(t')dt' 
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(6) 



where T signifies time-ordering. The correspond- 
ing eigenvalues of this operator are related to the 
quasienergies^ via Xj — exp[— iejT]. As the Floquet 
states provide a complete basis, the two-electron wave- 
function can be expanded as: 



!*(*)} =E a i exphM |^(t)). 



(7) 



P m n{t) = \(9(t)\m,n)\' 



(5) 



This expansion is directly analogous to the expansion of a 
wavefunction in energy eigenstates in the case of a time- 
independent Hamiltonian, and indeed in the adiabatic 
limit T — > oo the quasienergies evolve into eigenenergies, 
and the Floquet states to eigenstates. From Eq[7| it is 
clear that the long timescale behavior of the system is 
dictated essentially by just the quasienergies, since the 
Floquet states are periodic with period T, and also gen- 
erally have a weak time-dependence. Tunneling between 
two states is suppressed^ when their quasienergies ap- 
proach each other as an external parameter (such as the 
electric field) is varied. Thus as well as being directly ob- 
servable in the behavior of the electronic wavefunction, 
CDT may also be detected in the system's quasienergy 
spectrum by the presence of exact or avoided crossings. 
We shall show in the next section how these crossings 
in the exact quasienergies, obtained from the numerical 
diagonalization of the time-evolution operator © , corre- 
spond to the suppression of the dynamics of the electrons, 
and how analytic approximations to the quasienergies al- 
low the parameters at which this occur to be predicted 
with high accuracy. 



III. RESULTS 

In Fig^ we show the exact quasienergy spectra with 
parameters u = 2 and U — 16, for two extreme cases 
of array size: a two-site system (i.e. a double QD) 
and a large system of N = 16 sites. To assist the 
comparison between the two systems, the quasienergies 
have been classified into three categories according to 
how the electrons are distributed within the QD ar- 
ray in the corresponding Floquet state. To determine 




FIG. 1: (Color online) Quasienergy spectrum for (a) a two-site system, and (b) a 16-site system, for U = 16 and ui = 2. 
Symbols indicate the characteristic of the corresponding Floquet state: hollow diamonds (red) = doubly-occupied states, solid 
diamonds (black) = neighbor-states, points (blue) = states with wide separation between electrons. The widely-separated 
states (points) exhibit a series of miniband collapses in excellent agreement with EqlHJover the whole range of E. For weak 
driving the doubly-occupied states show a similar set of collapses but with half the period. The vertical dashed line marks the 
boundary between the weak and strong driving regimes. 



this we first class the basis states as i) states in which 
both electrons occupy the same QD ("doubly-occupied 
states"), ii) states in which the electrons occupy neigh- 
boring QDs ("neighbor-states") and Hi) states in which 



the electrons are more widely-separated. We then assign 
each Floquet state to the class onto which it projects the 
most, according to the magnitude of the overlap integrals 
|(^(0)|m,n}| 2 . 



Some similarities between the two spectra in FigQ] 
are immediately clear. In the two-site case a pair of 
quasienergies make a sequence of exact crossings as E 
is increased from zero up to U, while the remaining 
quasienergy stays isolated from them. As the two inter- 
crossing states evolve from the excited eigenstates of 
the static Hamiltonian, they consequently project mainly 
onto the doubly-occupied states. The remaining, iso- 
lated, Floquet state evolves from the system's ground 
state and therefore has the character of a neighbor-state. 



In the 16-site system we can see a similar, but richer, 
behavior. Instead of just two inter-crossing states in 
the weak-field regime, the set of doubly-occupied states 
form a miniband, similar to that observed for non- 
interacting electrons in RefsASi. In analogy to the cross- 
ings seen in the two-site system, this miniband exhibits a 
sequence of "band-collapses" at which the quasienergies 
become degenerate. These collapses occur at the same 
values of E at which the crossings occur in the two-site 
case. This behavior is seen more clearly in Fig|2i for 
a nine-site system. Here seven of the doubly-occupied 
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quasienergies make up a miniband and exhibit a set of 
band-collapses as E is increased, while the two remain- 
ing states, which physically correspond to the two edge- 
states, are almost degenerate and show little dependence 
on E. 

For non-interacting particles, Holthaus^ obtained an 
approximate expression for the quasienergy spectrum: 

6j = ^ J (E/lu) cos akj (8) 

where Jo is the zeroth Bessel function of the first kind, 
A is the width of the miniband, and kj are the permit- 
ted lattice momenta (assuming periodic boundary condi- 
tions). In the system we consider, the Coulomb interac- 
tion only operates if both electrons occupy the same QD. 
Accordingly, if the electrons are distant from each other 
they behave effectively as free particles, and in Fig^jD we 
can see that the Floquet states in which the two electrons 
are widely separated (plotted with points) form a mini- 
band of exactly the type given by Eq|S] Clearly these 
states cannot arise in the two-site system, as it is im- 
possible for the two electrons to be further apart than a 
single lattice spacing. 

Surprisingly, Holthaus' result for non-interacting par- 
ticles can also be applied to understand the miniband 
behavior of the doubly-occupied states, in which the 
Coulomb interaction is maximized. In this CclSG, cLS the 
miniband is composed of Floquet states that project onto 
doubly-occupied states, it is plausible to visualize the lo- 
calized two-electron state as tunneling from site to site 
as a single object. Since this object has charge 2e the pe- 
riod of the miniband is consequently halved, taking the 
form e,j — (A/2) J q (2E/uj) cosak r In Fig [2^ we plot the 
envelope of this miniband, and see that for low values of 
E the agreement with the exact results is excellent. As 
the value of E increases, however, the doubly-occupied 
Floquet states begin to mix and interact with the other 
states, and the miniband structure starts to break down. 
When E exceeds the size of the Hubbard gap the struc- 
ture is completely lost, and the system enters the strong- 
driving regime. 

To illustrate the physical significance of the miniband 
collapses, we show in Fig|2t> the minimum value of the 
probability P^t) attained by the system during ten pe- 
riods of the driving, having been initialized in the state 
1 5 , 5) - that is, with both electrons in the central QD. 
It can readily be seen that at the points of miniband 
collapse the system remains frozen in its initial state, 
thereby manifesting CDT, and producing sharp peaks 
in -P m i n centered on the points of quasienergy degen- 
eracy. This demonstrates that interpreting the two- 
electron state as a single composite particle, as suggested 
by the form of the Floquet spectrum, indeed provides the 
correct description of the system's dynamics. 



To complete the comparison of the two cases in Fig^ 
we now consider the strong-driving regime (E > U). 




FIG. 2: Weak-field regime for a nine-site system, U — 16 
and uj — 2. (a) Seven quasienergies form a miniband, 
while the two almost degenerate edge-states remain isolated 
from the rest. The solid lines are the approximate solution 
±(A/2) Jo(2E/u), where A is the width of the miniband. As 
E approaches U the miniband structure progressively breaks 
down, (b) Minimum value of Pss(t) for this system, initial- 
ized in state |5, 5). At the points of quasienergy degeneracy 
(miniband collapse) CDT occurs. 

Again, clear similarities between the two quasienergy 
spectra are evident. In the case of the two-site system, a 
pair of the quasienergies make a set of close approaches 
while the third oscillates weakly around zero. Detailed 
examination of the close approaches (RefiiS) reveals that 
they are in fact avoided crossings, and a perturbative 
treatment of this system performed in that work yields 
the approximate solution: 

e± = ±2J n (E/u), (9) 

where ncu — U. When this resonance condition be- 
tween U and the driving frequency is not satisfied, the 
quasienergies in this regime show little structure, and 
CDT is in general not produced. In Fig^i it is clear that 
a generalization of this phenomenon occurs, and a simi- 
lar set of close approaches between neighbor and doubly- 
occupied Floquet states occurs. In Fig|3^ we plot the 
quasienergies of the neighbor and doubly-occupied Flo- 
quet states for a nine-site array in this regime, and find 
that the locations of the quasienergy degeneracies are 
again given accurately by the perturbative result (Eq|5J|, 
namely at the roots of Js{E/uj). It can also be seen that 
the agreement between the perturbative approximation 
and the exact result improves as E increases. This is to 
be expected, as the "perturbation" used in the derivation 
of Eq[!5]consists of the tunneling component of the Hamil- 
tonian, and thus the result becomes exact in the limit of 
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FIG. 3: (Color online) Strong field regime for the nine-site 
system, (a) Crossings occur between Floquet states of doubly- 
occupied type (hollow symbols (red)) and those of neighbor 
type (solid symbols (black)). The lines are the perturbative 
solution ±2Js(E/ui). (b) Minimum value of -Pss(t) for this 
system, initialized in state |5, 5). CDT again occurs when the 
quasienergies cross. 

large E when the tunneling is negligible in comparison 
with the ac-driving. 

In Fig|3|3 we confirm that the structure seen in the 
quasienergy spectrum indeed corresponds to a modifica- 
tion of the system's dynamics. As before, the system is 
initialized in state |5, 5), and integrated over ten cycles of 
the driving field. The sharp peaks again visible in P mm 
are located at the points of quasienergy degeneracy, indi- 
cating that CDT occurs as expected. Thus in both weak 
and strong driving regimes, Floquet analysis can be used 
to predict the field parameters at which CDT will occur 
in a driven QD array with high accuracy. 
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FIG. 4: Time-dependent projection of the two-electron state 
onto the basis functions (see Eq[!)J with the system initialized 
in state |5, 5). Only the projections onto the doubly-occupied 
basis states are shown, as all the others are negligible in size. 
Field parameters: E = 2.35, u = 2. The differences between 
P44 and Pee, and between P33 and P77, are too small to be 
seen here. 

the wave-function of a nine-site system, initialized with 
both electrons in the central QD. The strength of the 
electric field has been set to E = 2.35, close to the point 
of the first miniband collapse at E ~ 2.4048, and thus the 
effective tunneling is heavily suppressed but non-zero. It 
can be seen that over the first hundred periods of the 
driving the initial state |5,5) evolves chiefly into a su- 
perposition of the states |4, 4) and |6, 6), by diffusing 
symmetrical!}^ to the neighboring QDs, while retaining 
the two-particle superposition. This is confirmed by the 
fact that throughout the time-evolution, the total projec- 
tion of the wavefunction onto the doubly-occupied states 
never takes a value of less than 0.955, indicating that the 
role of the other states is merely to act as short-lived 
intermediates in the diffusion process. 



We now examine the dynamics of the system in the 
weak-driving regime in more detail. If the system is ini- 
tialized in a state in which both electrons occupy the 
same QD, one would expect that the Coulomb repul- 
sion would act to repel the two electrons away from each 
other. We have already seen, however, that this does not 
occur, and that in the presence of the ac-field the two 
electrons in fact bind together and propagate as a single 
composite particl o 14 i 15 . Although surprising, the ability 
of an ac-field to convert a repulsive interaction into a 
time-averaged attraction forms the basis of an increas- 
ing number of applications, ranging from the well-known 
Paul trap^i to a recent proposal to bind atomic nuclei 
together— with laser fields. We show in FigQ] the ex- 
plicit time-dependence of the dominant components of 



Fig vividly demonstrates how the initially localized 
two-electron state on site 5 splits smoothly into two 
pulses moving to the left and right at an almost constant 
velocity. Although reminiscent of soliton dynamics, the 
lack of any dispersive effects to provide "refocusing" of 
the wavefront means that as the pulses traverse the sys- 
tem they do smear out to an extent. A result of this, 
for example, is that the peak in P33 does not rise to 
the same height as P44 in FigQ] due to blurring of the 
pulse-shape. Even so, the definition of the charge pulses 
is sufficiently crisp to suggest the possibility of using an 
array of QDs as a type of coherent electron turnstile 24 . 
In such a device the ac-field is used to control the rate 
at which the charge-pulses propagate, so that in a given 
interval of time a specified quantity of charge is trans- 
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FIG. 5: Time-dependent charge- density, |^(t)| 2 , for a nine- 
site lattice, with the initial condition of both particles on site 
5 (the central site). White / black signifies high / low values 
of density. Field parameters are E = 2.35, u> — 2. 

ferred along the array. The well-defined trajectories of 
the current pulses visible in Fig[S] indicate that the ve- 
locity of the pulse can be unambiguously measured, and 
in FigElwe give a plot of this quantity, obtained as the 
reciprocal of the time required for the initial QD to be- 
come empty. Since the velocity of the pulse is directly 
related to the magnitude of the inter-dot tunneling, the 
velocity shows the same Bessel function dependence as 
the quasienergies, and vanishes at field parameters such 
that Jq(2E/uj) = 0, at which the tunneling is completely 
suppressed. This clearly shows how the ac-field can be 
used as a control parameter to regulate the speed of prop- 
agation of the pulse by means of CDT. 



It is interesting to note that this arrangement generates 
an unusual current of entangled, correlated electrons. As 
the initial state propagates symmetrically to the left and 
right of the QD array, there is an equal probability of 
finding an individual electron on either side of the cen- 
tral QD, but due to the combined effect of the ac-field 
and the Coulomb interaction there is a very high chance 
that the other electron will also be occupying the same 
QD. Recently Saraga and Loss2^ have proposed using an 
array of three QDs to create currents of entangled elec- 
trons for use in quantum computing and communication 
applications. We note that their setup contains many 
similarities with the system we study here: a pair of elec- 
trons is injected into the central QD, which then emits 
electron pulses to its neighbors. The major difference is 
that in Ref— the current is composed of single electrons 



FIG. 6: Circles indicate the velocity of the charge pulse in the 
nine-site system as it propagates through the array from the 
initial state |5, 5). The solid line is vo\ Jo{2E/uj)\ where vo is 
a fitting parameter. 



(in a superposition of up and down spin states), whereas 
in our case the current is carried by a two-electron bound 
state. In this sense the driven QD array that we study 
can be regarded as an electronic analogue of an optical 
beam-splitter—, which divides the initial state into two 
spatially separated and entangled electron pairs moving 
at a velocity controlled by the ac-driving field. 



IV. CONCLUSIONS 

In summary, we have studied the dynamical behav- 
ior of two interacting electrons in an QD array under 
ac-driving. We have seen that the competition between 
the Coulomb energy and the driving field produces a rich 
phenomenology, which can be described well by means of 
Floquet theory. Throughout the weak field regime, the 
driving field and the Coulomb interaction combine in a 
counter-intuitive way to bind the two electrons together 
to form a composite particle resembling an exciton. The 
tunneling of this particle from site to site can then be 
controlled by manipulating the driving field, which per- 
mits either the complete freezing of the dynamics at the 
points of "miniband collapse", or the fine regulation of 
the current of entangled electron pairs. When the field 
is stronger than the Coulomb interaction the character 
of the Floquet states abruptly changes, and a different 
form of CDT occurs, arising from the close approaches of 
states of doubly-occupied type and neighbor- type. This 
generalizes the form of CDT observed previously in a two- 
site system^, and only occurs when the driving field is in 
resonance with the Coulomb energy (i.e. tilj — U) at the 
roots of J n {E/uS). The delicate control over the electron 
dynamics afforded by these effects has many potential 
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applications to the coherent manipulation of correlated 
electron states in mesoscopic systems. 

Although in this work we only consider one specific 
set of system parameters and one driving frequency, we 
have found that the results are of general validity as 
long as U > t (which is usually satisfied for QDs in the 
Coulomb blockade regime) and that the system is in the 
high-frequency regime 2 (u> > t). A recent experimental 
studj^ gave estimates of U ~ 3.7 meV and i ~ 0.5 meV 
for the double QD system studied there. Given that the 
QD separation was about 200 nm, the required driving 
field to see the CDT effects we report should be of the 
order of 200 GHz, with field-strengths of up to 10 kV/cm. 
Such systems also possess excellent coherence properties, 
with decoherence times of the order of nanoseconds re- 
ported recently^. We therefore believe that experimen- 
tal observation of these effects is feasible, and could be 



achieved, for example, by connecting the system to exter- 
nal leads and measuring the current passing through the 
array, or more directly, by measuring the charge occupa- 
tion of the individual QDs by means of quantum point 
contacts, as has recently been achieved in Refsiii 2 ^. This 
latter technique has the important advantages that the 
occupation of the QDs can be determined even if the 
inter-dot current is too small to measure by conventional 
means, and also should induce less decoherence than the 
transport form of measurement. 
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